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Abstract 

We study the growth of matter perturbations beyond the hnear level in cosmologies in which 
the dark energy has a variable equation of state. The non-linear corrections result in shifts in the 
positions of the maximum, minima and nodes of the spectrum within the range of Baryon Acoustic 
Oscillations. These can be used in order to distinguish theories with different late-time variability of 
the equation of state. 

PACS numbers: 95.36.-Fx, 95.35.-Fd, 98.80.Cq 



1 



1 Introduction 



The evolution of the matter power spectrum as a function of redshift depends on the underlying 
cosmological scenario. Its form reflects the evolution of the Universe since the time of matter-radiation 
equality. For given initial conditions, determined by the primordial spectrum (usually assumed to be 
scale invariant), the growth of matter perturbations can be used in order to constrain the possible 
cosmological scenaria through the comparison of the resulting spectrum with the observed large-scale 
structure. The most promiment feature of the matter spectrum is a series of peaks and valleys, 
characterized as Baryon Acoustic Oscillations (BAO). They originate in the period of recombination, 
and correspond to sound waves in the relativistic plasma of that epoch. 

The chacteristic length scale of BAO is around 100 Mpc. The exact form of the matter power 
spectrum at such scales is not easy to compute precisely, because of the failure of linear perturbation 
theory to describe reliably the growth of the corresponding fluctuations under gravitational collapse. 
At length scales below about 10 Mpc, the evolution is highly non-linear, so that only numerical 
N-body simulations can capture the dynamics of the formation of galaxies and clusters of galaxies. 
However, fluctuations with length scales of around 100 Mpc fall within the mildly non-linear regime, 
for which analytical methods have been developed. We focus on scales in the range 50-200 Mpc, 
within which BAO are visible. 

The various analytical methods [I]-[10] that have been developed in order to go beyond linear 
perturbation theory essentially amount to resummations of subsets of perturbative diagrams of ar- 
bitrarily high order, in a way analogous to the renormalization group (RG). (For a summary and 
comparison of the various methods, see [H].) We follow the approach of [5j, named time-RG or 
TRG, which uses time as the flow parameter that describes the evolution of physical quantities, such 
as the spectrum of perturbations. The method has been applied to ACDM and quintessence cosmolo- 
gies [5], allowing for a possible coupling of dark energy to dark matter [6^, as well as models with 
massive neutrinos [7]. 

The purpose of the present work is to apply the formalism to the case of a variable equation of 
state of the dark energy component. The typical example of such a scenario links the dark energy to 
an evolving scalar field with non-zero potential [12] ■ We use a simpler parametrization of the dark 
energy sector, by assuming that its energy-momentum tensor has the usual perfect-fluid form but 
with an equation of state that is a function of redshift: p/p = w^z). We assume the form of wlz) 
employed in [15] : 

aw{a) atrans 

w{a) = — , (1) 

a + atrans « + At rans 

where a — + z) is the scale factor, and w(a) = wo + {1 ~ a)wa- The value atrans corresponds to 
the "transition epoch" , during which the dark energy ceases to have the form of a pure cosmological 
constant with ui = — 1 and starts to evolve. At small redshifts eq. ((TJ reproduces the commonly used 
linear form of w(a) [14] 

w{a) = 11)0 + {I — a)wa- (2) 
Following [13], we describe the various models through the present-day value of w, wo = w{z = 0), 
its first derivative, w' = dw/dz\z^o, and atrans. For atrans <C 1, we have wo — wo and w' ~ Wa- A 
detailed discussion of this parametrization is given in appendix C of [13] . 

For our study we use values of atrans such that the equation of state is very close to that of a 
cosmological constant during the time of recombination. In this way, we do not account for potentially 
interesting models with significant amounts of evolving dark energy at early times. Instead, the 
emphasis of our study is on the infiuence of the late-time evolution on the matter spectrum. Assuming 
that the dark energy density is constant during recombination, with a negligible contribution to the 
total energy density, provides the most convenient framework in order to isolate the late-time effects. 



2 Formalism 

We concentrate on the evolution of fiuctuations at subhorizon scales. We include contributions 
from all massive matter (baryonic and dark) in the total matter density p(t,x). The fundamental 
equations in the TRG approach are the continuity, Euler and Poisson equations. For an expanding 
background, they can be obtained starting from the Einstein equations and the conservation of the 
energy-momentum tensor. This calculation has been performed in [B] in the case that the dark energy 
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arises through the potential of an evolving scalar field (j). We summarize the results, in order to deduce 
the relevant equations for the parametrization through a variable equation of state w{z). The details 
are given in the appendix A of ref. '6'. 

For the metric, we consider an ansatz of the form 

ds^ = a{T) [(1 + 2$(r, x)) dr^ - (1 - 2$(r, x)) dx dx] . (3) 

We assume that the Newtonian potential $ is weak, $ <C 1, and that the field can be decomposed 
as (;A(r, x) — (j>{T) + 5ij){T, x), with 5cf>/cj) <^1. In general, = 0(1) in units of the reduced Planck mass 
M = (SttG)^^''^. In the absence of a coupling between dark matter and dark energy, the magnitude 
of the fluctuations of (f> is not expected to exceed that of the gravitational field $. Finally, the density 
can be decomposed as p{t,x) = p(r) + 5p(T,x). We allow for significant density fluctuations, even 
though our analysis is not applicable when they are much larger than the background density. 

For subhorizon perturbations with momenta k ^ T-L = d/a, the linear analysis predicts \5v\ ~ 
{k/T-L)^ ~ {'H/k){Sp/p). A consistent expansion scheme can be obtained if we assume that $ <C 
\6v\ <SC 1. Including the density perturbations, our assumptions can be summarized in the hierarchy 
of scales: $ <^ \Sv\ <^ Sp/p < 1. We assume that the fluctuations of the scalar field 5(j)/(j) can be as 
large as the gravitational potential However, we shall see in the following that, if there is no direct 
coupling between dark matter and dark energy, they are actually much smaller. As we are dealing 
with subhorizon perturbations, it is consistent to make the additional assumption that the spatial 
derivatives of 5<j) dominate over their time derivatives. The predictions of the linear analysis allow 
us to make a more quantitative statement. We assume that a spatial derivative acting on 5(j) or 5v 
increases the position of that quantity in the hierarchy by one level. In this sense V"l> is comparable 
to 5v, while is comparable to p. 

The equation of motion for the scalar field has been derived in ref. [6]. Assuming the hierarchy 
we described above, it can be split into an equation for the homogeneous part: 

| + 2^^ + a^^ = 0, (4) 
a d<p 

and one for the perturbation: 

S;i> + 2^50 - V^50 + a'^W^H - m + = 0. (5) 

a d4>-' dtp 

The terms with time derivatives are subdominant according to our assumptions. The scalar field (j) 
takes values of order M. For (j) to evolve at cosmological times, it must have a mass term d^U{<f)) /dcf)^ = 
OiU^). We also have U{^), p = 0{'H'^M^). It is natural to expect dU{(j>)/d(t) = 0{'H'^M). For 
subhorizon perturbations with characteristic momenta k ^ "H, the inspection of eq. (O demonstrates 
that 5(j>/(j) is suppressed by (H/k)^ relative to In the absence of a direct coupling between dark 
matter and the field 0, the field fiuctuation 50/0 is negligible within the hierarchy we have assumed. 

We can conclude that the fluctuations in the dark-energy sector do not affect the growth of matter 
fluctuations at scales much smaller than the horizon distance. If the dark energy does not couple 
directly to matter, it does not cluster sufficiently for its fluctuations to play a role in structure for- 
mation. This conclusion is consistent with the approximations underlying the essentially Newtonian 
framework that we employ for the study of the spectrum. 

The equation of motion for the gravitational potential $ can be obtained from the flrst Einstein 
equation. The leading terms give 



a J 3M2 



(6) 



while our assumptions reduce the equation for the perturbation to the Poisson equation for the 
gravitational field $: 



V^<J> = ^a^Sp. (7) 
We now turn to equations derived from the conservation of the total energy-momentum tensor 
r^ff ~ 0- F'^'^ M = 0, the leading terms give 



-p + snp = ~l(^<i> + 2n$ + a'^yo, (8) 
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where we have employed eq. Q. The equation for the perturbation is more complicated. It can be 
simplified considerably through our assumptions about the hierarchy of the various fields. We obtain 



6p + m5p + V[{p + 5p)5v\ = -i^V'5(j)~Q, (9) 



where we have made use of our earlier conclusion that 5</>/^ is suppressed by (H/k)'^ relative to $, 
and therefore by {'H/k)*' relative to 5p/p. For = i = 1, 2, 3 we obtain the generalization of the Euler 
equation for this system. After eliminating higher-order terms and employing eqs. Q, ((8]), (O we 
find 

5ij+'H5v+ (SvV\Sv = ~V^. (10) 



We have seen that, if the dark energy has no direct coupling to matter, it influences the matter 
spectrum only through its effect on the expansion rate. We can now rewrite the relevant equations 
using the variable equation of state w(a). We replace the equation of motion of cj> by the continuity 
equation for the dark-energy density in the form 

PDE + 'i{l + w{a))HpDE ^Q- (11) 

For a specific quintessence model, the determination of the function w(a) requires the integration of 
eq. (|4|. In our study, however, we shall assume that w(a) is parametrized by eq. The Friedmann 
equation ((BJ takes the form 

H^=(^)^^a^(p-.P..). (12) 

In summary, the evolution of the background is determined by eqs. ([S]), Hll[) . (|12|l . while the growth 
of the matter fluctuations is described by eqs. ((T)), (O, (HOJ. 

The evolution equations are expressed in their most useful form in terms of the density contrast 
(5 = 5p/p < 1 and S(k, r) = V ■ 5ti(k, r). For the Fourier transformed quantities, we obtain 

<j(k, r) + 6'(k, r) j d^kid^k2<5i5(k-ki -k2)a(ki,k2)6l(ki,r)5(k2,r) = 0, (13) 

with 3(ki,k2) — ki ■ (ki +k2)/fcf, and 

e(k,r)+W6l(k,r) + |?^2nM<5(k,r) + j d^ki d^k2 (k-ki -k2) /3(ki, k2) 6l(ki, r) e(k2, r) = 0, (14) 

with /3(ki, k2) = (ki + k2)^ki • k2/(2fc?fei) and Q,m = pa^/{m^). 
We define the doublet 

/ V5i(k,»?) \ 




where rj — In a(r). This allows us to bring eqs. (|13p . (|14p in the form [TJ|?1[5] 



(15) 



a^iy3a(k,77) + Q,ab'^b{^,'n) = e^-yabci^, -ki , -k2 ) (/Pft (ki , r;) 1^3^ (k2 , r?) . (16) 

The indices a, h, c take values 1, 2. Repeated momenta are integrated over, while repeated indices are 
summed over. The functions 7, that determine effective vertices, are given in [4] [5]. The non-zero 
components are 

7i2i(k,ki,k2) = ^^^^^^<5D(k + ki +k2) = 7ii2(k,k2,ki) 

2 (17) 

7222(k,ki,k2) = ^(ki,k2) 5i3(k + ki + k2). 
The fi-matrix, which determines the linear evolution, is 

«(^) = I 3 ^, h (18) 

^-^M 2 + - 
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where a prime denotes a derivative with respect to i]. 

The next step is to derive evolution equations for the power spectra. The spectrum and bispectrum 
are defined as 



{</?a(k, »7)(/56(q, n)) =<5D(k + q)Pai,(k, i]) 
(VPa(k, r?)<(5b(q, ri)Vc{p, v)) ^^(k + q + p)Ba6c(k, q, p, 77) 



(19) 



Their evolution with time can be obtained through differentiation with respect to rj and use of eq. 
(|16p . The essential approximation that we have to make in order to obtain a closed system is to 
neglect the effect of the trispectrum on the evolution of the bispectrum. This gives [H [5] 



dr,Babc{K -q, q - k) = -QadBdbcCk., -q, q - k) - ^ItdBadcCk., -q, q - k) - ^lcdBabd{K -q, q - k) 



The procedure of truncating the system of equations is commonly employed in the applications of 
the Wilsonian RG to field theory or statistical physics. (For a review, see [IS].) The accuracy of the 
calculation can be determined either by enlarging the truncated system (by including the trispectrum, 
for example) and examining the stability of the results, or by comparing with alternative methods. 
The second approach is often followed, because enlarging the truncation can increase the complexity 
of the calculation considerably. In the case of the TRG, the agreement with results from N-body 
simulations for ACDM [5], quintessence and coupled quintessence [6] is at the sub-percent level for 
the BAO range. We expect agreeement with alternative methods at the same level for the case of 
a variable equation of state we are studying. We shall see in the following that the features of the 
spectrum differ at the sub-percent level for various forms of the equation of state. We emphasize, 
however, that in this work we make comparisons between results obtained always with our method 
(e.g. for various forms of 'w{z), at the linear or non-linear level), without reference to other approaches. 
As a result, the bigger systematic effects, the ones that generate the discrepancy between different 
methods, are expected to cancel out. In particular, for the parametrization of eq. ([T}, we expect that 
the comparison we perform between models characterized by different values of wo and w' can lead 
to the reliable identification of the relative variation of features in the spectrum, such as the location 
of maxima, minima or nodes. 

3 Results and discussion 

We consider models for which ^1(2) is very close to —1 at high redshfit. We also assume a common 
normalization at high redshift for the matter spectra of all the models that we study. At the linear 
level, the location of the extrema of the spectrum at low redshift is expected to be independent of the 
form of w{z). The reason is that, at the linear level, the evolution of the matter spectrum amounts 
to multiplication by an overall redshift-dependent factor (the linear growth rate). At the non-linear 
level, however, the shape of the matter spectrum is modified through mode coupling during the 
growth of the fiuctuations. The main purpose of our study is to determine the magnitude of the shift 
of the peak location arising from the non-linear corrections for the various forms of w{z). 

The full system of eqs. (|20|l . (|2ip can be solved in a way analogous to that described in Appendix 
B of [5]. We set the initial conditions for the integration of the evolution equations for the spectra at 
a redshift z = 40. At such early times the evolution is linear to a very good approximation. Moreover, 
our assumptions about the form of w{z) imply that there is no appreciable amount of dark energy at 
such early times. As a result, the early evolution of the spectrum is identical to that in the ACDM 
case. We compute it by making use of the numerical code CAMB [16]. We assume that the primordial 
spectrum is scale invariant with spectral index n — 0.96. We use the Cosmic Microwave Background 



dnPabiKv) 



I d^q[^acdO^: -q, q - k)B6cd(k, -q, q - k) 




(20) 




(21) 
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Figure 1: Linear and non-linear spectra at z = 
0, for Wo = —0.8, w' — —0.7, normalized with 
respect to a smooth spectrum. 



Figure 2: The fractional shift of the maximum, 
minima and nodes of the non-linear spectrum, as 
a function of redshift, for wq = —0.8, w' = —0.7. 



(CMB) normalization for all the models we study: = 2.46 x 10"'^ at feo = 0.002/Mpc. We take 
the present-day Hubble parameter to be Hq = 70 km/sec/Mpc, the present-day total matter density 
Hm — 0.27 and the baryonic density Q,t — 0.046. Under the assumption that the dark energy is 
negligible before z = 40, all the models that we have studied have the same early Universe history. It 
is, therefore, consistent to assume that they have the same power spectrum at z = 40. The differences 
appear at low redshifts at which the equation of state w{z) deviates from a pure cosmological constant 
with w = —1. We take atrans ~ 1/(1-1- ztrans), with ztrana = 10, for the "transistion epoch", during 
which the dark energy starts deviating from a pure cosmological constant. The various models predict 
different present-day spectra, age of the Universe, amplitude of matter fluctuations (crs), and location 
of BAO peaks. 

In flg. [T]we depict the linear (solid line) and non- linear (dashed line) matter power spectra at 
z — for a model with wq = —0.8 and w' = —0.7. Following the commong practice, we have 
normalized both spectra with respect to a smooth spectrum. (The corresponding transfer function is 
given by eq. (17) of ref .17 .) We have used a smooth function such that the non-linear spectrum is 
oscillatory around an almost constant value in the fc-range of interest. On the same figure we have 
indicated certain characteristic points of the spectrum: two minima (I and 5), a maximum (3) and 
three nodes (2,4 and 6). In fig. [2] we depict the evolution of these points as a function of redshift. 
At the linear level, no evolution is expected. For this reason, we plot the fractional difference in the 
location of these points from their constant values at the linear level. We observe that for large z 
the deviation is negligible, because the spectrum is linear to a good approximation. The deviations 
become significant below a redshift approximately equal to 3. We have checked that the choice of the 
smooth spectrum used for normalization does not infiuence appreciably the shift of the characteristic 
points. 

The change in the shape of the spectrum arising through the non-linear corrections can be used in 
order to distinguish different models. The ACDM maximum /cs^acdm is shifted by 6.4 x lO'^/i/Mpc 
when non- linear corrections are taken into account. The shift is different for models with a variable 
equation of state. In fig. Owe depict the location of the maximum of the non-linear spectrum as 
a function of wo and w' at a redshift z = 0.366. We plot the fractional deviation of the maximum 
from its position in ACDM (corresponding to the point wo — ~1, w' — 0). We display several points 
that correspond to distinct theories, as well as a smooth fit in order to guide the eye. There are no 
points beyond the line starting at (—0.6,0.5) and ending at ( — 1.5, 1.5). In that range of the {wo,w') 
plane, and for our choice of atrana — 0.091, the values of the Hubble rate H ai z — 40 differ by 
more than one per mille between the various models. This implies that the effective equation of state 
w{z) deviates from that of a pure cosmological constant already at the initial time at which we have 
started the evolution of the spectrum. As we wish to disentangle such early-time effects from pure 
late-time ones, we do not investigate this part of the parameter range. It is apparent from fig. [3] that 
the late-time non-linear corrections to the spectrum generate differences at the per mille level for the 
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Figure 3: The fractional shift of the first maximum from its location for ACDM, as a function of wq and 
w', at a redshift z = 0.366. 

location of the maximum in various dark energy models. 

Similar behavior is observed for the location of the minima and nodes of the spectrum. The ratio 
of the position of the first minimum k\ to that of the first maximum fca , as a function of wq and w' at 
a redshift z = 0.366, results in a surface with shape very similar to that in fig. |3] but with opposite 
tilt. The reason is that the first minimum is located in the part of the spectrum that is described 
accurately at the linear level, so that k\ is essentially the same for all the theories we consider. In 
table [T] we give the values of fci/fca for ACDM and four theories that correspond to the corner points 
in fig. 131 The ratio of the position of the second minimum to fca displays enhanced variation, 
because the non- linear corrections increase the value of k-j, but reduce the value of fcs. In table [1] we 
give the values of ks/ks for the same five theories. 

There are certain degeneracies in the form of the non-linear spectrum, which are not immediately 
apparent in fig. [3] In order to display them we depict in fig. |4]the same data points as a function 
of the total time t4o between z — 40 and z = 0, which is equal to the age of the Universe in a good 
approximation. We normalize this time with respect to its value for ACDM. It is apparent that the 
shift of the maximum has an almost linear dependence on t4o. Models with distinct values of wo and 
w' , but the same value of t4o, predict a similar shift of the maximum. This indicates that, if Hm is 



Table 1: The location of mimima relative to the maximum for various theories. 
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Figure 4: The fractional shift of the first maxi- 
mum from its location for ACDM, as a function 
of the age of the Universe. 



Figure 5: Angular distance as a function of the 
age of the Universe for the class of models we 
consider. 



kept fixed in the various models, as we have assumed throughout this work, the growth of non-linear 
perturbations depends mainly on the available time. An equivalent measure is the value of erg, as 
computed from the linear spectrum for each model. For fixed ^Im, the value of as increases with the 
available time. This results in an almost linear dependence of the shift of the max;imum on erg. 

The various models that result in a shift of the maximum very close to that for ACDM lie roughly 
on a line that starts at (—0.6,-1.5) and ends at (—1.5,1.3) on the {wq,w') plane. The same also 
happens for the shift of the second minimum, whose deviation from the value for ACDM can be as 
big as 1%. The line between (—0.6, —1.5) and (—1.5, 1.3) lies within the best-fit range obtained from 
a combination of the WMAP 5-year data, supernova data and the linear BAO specturm (fig. Dl 
of ref. [13]). The reason can be traced in the strong correlation between the angular distance of a 
high-redshift source and the corresponding time. In fig. [5] we depict the angular distance of a source 
at z = 40 as a function of the total time t^o between z — 40 and z — 0, for a large number of models 
with various values of lOo and w' . The same relation also holds for higher redshifts, as all the models 
have the same evolution before z — 40. It is apparent that for the class of models that we consider 
the time determines almost uniquely the angular distance for high-redshift sources. As the time also 
determines the growth of the non-linear effects in the matter spectrum, the appearance of the same 
degeneracy in the CMB and matter spectra is natural. 

In conclusion, we have presented a calculation that demonstrates how the matter spectrum can be 
used in order to differentiate between cosmological models with a late-time deviation from ACDM. 
Features of the spectrum such as the location of extrema and nodes do not remain constant when non- 
linear corrections to the evolution are taken in to account. The differences between various models 
are small (at the sub-percent level). However, they are calculable and can be used for comparisons 
with the spectrum deduced from observations. The most efficient method would make simultaneous 
use of the full spectrum, as the positions of extrema and nodes vary differently. For example, non- 
linear corrections shift the location of the first maximum towards larger values, while they decrease 
the value of the second minimum of the spectrum. A simultaneous fit to the locations of maxima, 
minima and nodes of the spectrum can provide a useful method for differentiating between theories 
with late-time variation of the equation of state. 
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